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The Mc@Nlo and MePs@Nlo methods, as implemenated in the Monte-Carlo event generator 
framework Sherpa, are used to estimate the perturbative and non-perturbative uncertainties in 
various processes such as dijet production and the production of a W boson in association with 
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1. Introduction 

Being largely stimulated by the need for higher precission of theoretical predictions in both 
Standard Model analyses and new physics searches at the LHC, the simulation of higher-order 
QCD corrections in Monte Carlo event generators has seen vast improvements in recent years. To 
this end, two lines of development have been followed. In the MePs approach [1, 2, 3, 4, 5, 6] 
higher-order tree-level matrix elements of successive final state parton multiplicity are merged into 
an inclusive sample, offering both leading-order accuracy for the production of hard partons and re- 
taining the overall resummation of scale hierarchies through the parton shower at the same time. On 
the other hand NloPS approaches, introduced as either Mc@NLO [7] or POWHEG [8, 9], work on a 
single parton multiplicity elevating its accuracy to next-to-leading order. Both methods have been 
shown to be automatable [10, 11] within the Sherpa event generator framework [12]. Thereafter, it 
was sought to recombine both lines of development. In a first step, called the MEnloPS prescrip- 
tion, the NloPS and MePs methods have been combined using the NloPs' NLO accuracy for the 
inclusive process supplementing it with higher-order tree-level matrix elements in an MePs fash- 
ion [13]. In second step multiple NloPS processes of successive parton multiplicity are combined, 
elevating the accuracy of the MePs method to next-to-leading order, dubbed MePs@Nlo [14, 15]. 
In the following both the NloPS and MePs@Nlo methods are summarised. Particular emphasis is 
put on both methods' major accompishments with respect to standard leading order computations: 
its increased theoretical accuracy expressed through reduced perturbative uncertainties. 

2. NloPS matching 

Following the notation of [11] a general NLO+PS matching can be cast in the form of the 
following master formula 

(O) = j^ B B< A >(*|,) A (A) (to, 4) 0(<P B ) + Z j* °' 1 } A (A) (*, 4) <?(**) 
+ J d<D« H(<J>«) 0(® R ) . 

(2.1) 

Therein, the NLO-weighted normalisation of the resummed events is defined as 

B (A) (O b )=B(<D b ) + V(<D b ) + I (A) (<D b )+£ ( d4»! [D^ A) 0(^_ ? )_ D ( S) j (d**,^). (2.2) 

t = t(<&\) identifies the infrared limits of the additional parton's phase space and serves as an 
ordering variable of the parton shower resummation. The resummation kernels are then defined by 
the auxiliary set of subtraction kernels D (A) , ensuring the correct behaviour in both the soft and the 
collinear limit of the emission of an extra parton, exhibiting full colour and spin correctness in the 
respective limits. They imply the modified Sudakov form factor 
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An upper scale \Xq limits the region of resummation, i.e. the exponent of the Sudakov form 
factor vanishes at t = \Iq. This scale has been made accessible for the first time in the implemen- 
tation of [11] and can thus be used to study the uncertaity related to its arbitrariness. The finite 
remainder of the real emission cross section is then embedded in the so-called hard events denned 
through 

H(<J>*) = R(<*>«) - YPi K \^R) ■ ( 2 - 4 ) 

i 

Fig. 1 now shows an evaluation of the resummation scale uncertainty in various Mc@NLO imple- 
mentations for pp — > W + «jets [19] and contrasts it with the renormalisation and factorisation scale 
uncertainties in a standard fixed-order next-to-leading order calculation. Fig. 2 details all sources 



Rapidity Distance of Leading Jets Azimuthal Distance of Leading Jets 




-4-3-2-101234 o 0.5 1 1.5 2 2.5 3 

Ay(First Jet, Second Jet) A^fFirst Jet, Second Jet) 



Figure 1: Rapidity (left) and azimuthal (right) separation of the two leading jet pp — >> 2jets compared to 
ATLAS data [16]. 
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Figure 2: Left: Inclusive jet cross section in pp — »> 2jets compared to ATLAS data [17]. Right: 3-jet over 
2-jet ratio in dependence on the scalar transverse momentum sum of all jets in pp —>> 2jets in comparison 
to CMS [18]. 
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of perturbative }ip, }Xq) as well as non-perturbative uncertainties due to the multiple interac- 
tion model in an Mc@NLO implementation of inclusive and dijet production [20]. In all cases, the 
perturbative uncertainties for observables described at NLO accuracy are greatly reduced while the 
parton shower resummation provides the correct description when large hierarchies of scales in t 
are present. At the same time, there are observables/regions where the uncertainty on the modelling 
of the soft structure of the event is non-negligible. 



3. MePs@Nlo merging 

The NloPs matched calculations detailed in the previous section can now be used as input to 
extend the CKKW-type to next-to-leading order [15, 14]. The master formula for its construction 
reads as follows 

r 

p _ T r D (A) 

(O) = J d<D„B< A ) A ( ^(t c ,^ Q )O n + J d<J>! -^-A ( n A) (t n+i ,Li 2 Q )e(Q cut -Qn + i)O n+i 

tc " 

+ J d<D ;i+1 ni A) A ( ^ s \t n+1 ,^)&(Q cut -Q n+l ) O n+l 
( ^ \ 

+ J d$ n+1 5^ / d<j> lKn y^\t n+ul iim Qn+l -Q cux ) (3-D 



hi- I 



Y D (A) 

x A ) ( ! A) 1 (f c ,f„ +1 )O n+1 + / d*i ^-Af^{t n+2 ,t n+ i)O n+2 

J a n+l 

+ J d<J> ;i+2 H^A^(f„ +2 ,?„ +1 )Ai PS) (?„ +1 ,^)0(a i+ i -Gcut) O n+2 + ... , 

Therein an Mc@NLO description of an n parton multiplicity is restricted to have its emission 
produced at a jet measure Q smaller than Q cut . The region with Q > Q cut is then filled with an 
Mc@Nlo for the n + l parton process. To restore the correct resummation with respect to the n 
parton process to at least parton shower accuracy its Sudakov form factor A ; ( j PS) is inserted. The 
overlap with similar terms in is removed with the term in the braces on third line. A multijet 
merged description is then achieved by iteration eq. 3.1. 

Again, the calculation benefits from the decreased theoretical uncertainty of its Mc@NLO 
input processes. Figs. 3 and 4 exemplify this feature for the process pp — > W + jets compared to 
ATLAS data. For this calculation the processes with 0, 1 and 2 additional jets are described at next- 
to-leading order while 3 and 4 addiotional jets have been merged on top of that at leading order 
accuracy. These different levels of accuracy can be directly seen in the respective uncertainties. 
Further, they are contrasted with a MEnloPS [13] prediction using an Mc@NLO input only for the 
pp — » W process and merging only leading order prediction for 1, 2, 3 and 4 additional jets on top. 
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Figure 3: Cross section as a function of the inclusive jet multiplicity in pp —> W +jets events compared to 
ATLAS data [16]. 
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Figure 4: Differential cross section as a function of the transverse momentum of the first (left) and second 
(right) jet in pp —> W+ > l,2,3jets events compared to ATLAS data [16]. 
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